Run any of the following code
- POtest.m         % require CVX to produce a reference solution
- POtest_no_cvx.m  % if you do not have CVX installed
